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SUMMARY 

Understanding the relationship between taxonomic and 
morphological changes is important in identifying the rea- 
sons for accelerated morphological diversification early in 
the history of animal phyla. Here, a simple general model 
describing the joint dynamics of taxonomic diversity and 
morphological disparity is presented and applied to the 
data on the diversification of blastozoans. I show that 
the observed patterns of deceleration in clade diversifica- 
tion can be explicable in terms of the geometric struc- 
ture of the morphospace and the effects of extinction and 
speciation on morphological disparity without invoking 
major declines in the size of morphological transitions 
or taxonomic turnover rates. The model allows testing 
of hypotheses about patterns of diversification and esti- 
mation of rates of morphological evolution. In the case 
of blastozoans, I find no evidence that major changes in 
evolutionary rates and mechanisms are responsible for the 
deceleration of morphological diversification seen during 
the period of this clade's expansion. At the same time, 
there is evidence for a moderate decline in overall rates 
of morphological diversification concordant with a major 
change (from positive to negative values) in the clade's 
growth rate. 

1. INTRODUCTION 

Rapid morphological diversification early in a clade's his- 
tory, at relatively low taxonomic diversity, with an appar- 
ent slowdown afterwards represents a commonly observed 
pattern of radiation of animal life. Among the best known 
examples are Paleozoic blastozoans (Foote 1992; Wag- 
ner 1995a), bryozoans (Anstey & Pachut 1995) and gas- 
tropods (Wagner 1995b), Paleozoic and Mesozoic crinoids 
(Foote 1994, 1995, 1996b), Cambrian marine arthropods 
(Briggs et al. 1992; Wills et al. 1994) and Ordovician 



trilobites (Miller & Foote 1996). The pattern of rapid 
initial increase in morphological disparity which remains 
unsurpassed during the history of the clade afterwards 
has often been interpreted as an evidence for major secu- 
lar changes in evolutionary rates and mechanisms (Valen- 
tine 1969, 1980). Different explanations for these secular 
changes have been proposed. It has been argued that 
ecological opportunities were greater in the early history 
of many clades, that genetic and developmental systems 
were less canalized early on, and that the nature of adap- 
tation on a "rugged" adaptive landscape results in a slow- 
down of the rate of adaptation (e.g., Erwin 1994; Erwin et 
al. 1987; McShea 1993; Valentine 1969, 1980; Valentine et 
al. 1994; Kauffman 1993). Each of these factors can po- 
tentially cause a reduction in the probability and/or size 
of morphological changes with time which will translate 
into a decline in the rate of clade's diversification. 

In spite of the well-recognized potential importance of 
the pattern and extensive discussions of its generality and 
possible explanations, there have been only few attempts 
to use formal mathematical models to identify a minimum 
set of factors sufficient for explaining the pattern and 
to test hypotheses about underlying mechanisms. Ex- 
isting time-homogeneous models have predicted a linear 
increase in morphological variance (Slatkin 1981; Foote 
1991, 1996a; Valentine et al. 1994) reinforcing the be- 
lief that something has to change significantly in time 
to result in the observed patterns. However, some dif- 
ficulties arise when one tries to apply these models to 
data. One reason is that while the models consider sin- 
gle traits that vary continuously and whose evolution is 
unconstrained, the empirical studies of morphological dis- 
parity are typically based on a large number of discrete 
characters which are always subject to some morphologi- 
cal constraints (geometric, structural, or functional limits 
on possible trait values) . Another reason is that previous 
modeling frameworks did not include some of the fac- 
tors that can significantly affect the dynamics of morpho- 
logical disparity (such as subclade extinction or origina- 
tion events that do not result in large differences between 
sister-species). Here, I extend the previous work by con- 
structing a more detailed model of clade diversification 
specifically designed for treating discrete characters and 
by applying it to the data on the diversification of blas- 
tozoans. 

2. MODEL 

I consider the evolution of a monophyletic clade driven 
by extinction, speciation and anagenetic changes. Let 
us assume that each lineage in the clade is character- 
ized by L binary morphological traits. A lineage a can 
be described by a sequence of 0's and l's of length L: 
l a = {If , 1%, . . . , If) where 1% = or 1 (Jfe = 1, . . . , L). 
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The morphological space is mathematically equivalent to 
a binary hypercube. In discussing the clade's diversifica- 
tion, it is useful to visualize each lineage as a point on 
a vertex of the morphological hypercube. Accordingly, a 
clade will be a cloud of points. Speciation, extinction and 
anagenesis change the size, location and structure of this 
cloud. Let us define a morphological "distance" between 
lineages a and (3 as the number of traits at which the 
lineages are different: 

d a0 = J2(l?-lf) 2 - (1) 

I will be interested in the joint dynamics of the number of 
lineages in the clade, N, and the average pairwise distance 
within the clade, D = J2i<j (P j /(N(N - l)/2). Distance 
D is a measure of morphological disparity characterizing 
the spread of the clade in the morphological space. 

The dynamics of the clade size and the morphologi- 
cal disparity have been the focus of the previous work 
(e.g. Slatkin 1981; Foote 1991, 1996a; Valentine et al. 
1994). Here in addition to N and D, I will consider the 
average distance of the members of the clade from its 
species-founder, d. Let e denote the species-founder of 
the clade. Then d — J2i d el /N . Distance d characterizes 
the extent of the evolution of the clade from its ancestral 
state. Below this measure will prove to be very infor- 
mative and convenient to use in analyzing real data. I 
will model clade evolution as a random walk on the mor- 
phological hypercube with births and deaths. That is I 
will assume that origination and extinction events as well 
as morphological changes can be considered as random 
and independent of the morphology (cf. Raup & Gould 
1974; Slatkin 1981; McKinney 1990; Foote 1991, 1996a; 
Valentine et al 1994). This represents a null hypothesis 
which must be rejected before introducing additional fac- 
tors to explain the observed patterns of taxonomic and 
morphological diversification. 

It is convenient to formulate the model in discrete time. 
I consider two types of morphological changes: anage- 
netic and cladogenetic. Anagenetic changes are modeled 
by assuming that during a unit time interval each trait in 
a lineage may evolve to an alternative state with a small 
probability fix. I assume that there are two types of origi- 
nation events having probabilities a± and 02 , respectively. 
The origination events of the first type do not result in any 
immediate differences between two (or more) new lineages 
the old lineage has split into. This might be the case when 
different large parts of a subdivided populations become 
completely isolated by geographic (as in the vicariance 
speciation scenario, e.g. Lynch 1989) or reproductive (as 
in the parapatric speciation scenario, e.g. Gavrilets 1999; 
Gavrilets et al. 1998) factors. The origination events of 
the second type are accompanied by (significant) morpho- 



logical changes (Eldredge & Gould 1972). This might be 
the case when speciation takes place in a small (periph- 
eral) population which has undertaken significant mor- 
phological evolution before emerging as a new species (as 
in the peripatric speciation scenario, e.g. Mayr 1963). 
During such speciation events each trait in a new lineage 
can evolve to an alternative state with probability fi2- I 
assume that there are two types of extinction events. A 
lineage can become extinct individually (with probability 
S\) or as a member of a "T-subclade" simultaneously with 
all other members (with probability 62). Following Der- 
rida & Peliti (1991), I say that two lineages belong to the 
same T-subclade if their last common ancestor existed T 
years ago. Extinction of a subclade might happen when 
some traits that are shared by the members of the sub- 
clade "promote" extinction (McKinney 1997) (say, after 
a change in the environment). I will assume that all rates 
defined above are small (/Zj, <7j, Si << l,i = 1,2). 

The changes in N, D and d between subsequent time 
intervals are described by a system of difference equations 
(see Appendix for details) 

AN = (a- S)N, (2a) 
AD = - Ufi +^-+ D + 2fiL, (2b) 

Ad = -2(i(d- -). (2c) 

Here a = a\ + (j% and S = 8% + #2 are the overall rates 
of speciation and extinction (per unit time interval per 
lineage), and 4> is the proportion of the clade represented 
by a T-subclade that goes extinct. In general, these rates 
can change in time and/or with the clade size. With fos- 
sil data, in practice, it will generally not be possible to 
distinguish the two types of speciation events, and it will 
be very difficult to distinguish anagenetic from cladoge- 
netic morphological change, but the overall rates of ex- 
tinction and origination can be estimated. Parameter 
fx = fii +CT2M2 is the overall probability of a morphological 
change (per trait per unit time interval per lineage) which 
incorporates both morphological and taxomonic rates of 
evolution. Below I describe a simple method for estimat- 
ing fi from fossil data. 

The system of difference equations (||) describing the 
joint dynamics of N, d and D, can be easily solved nu- 
merically. Below I consider several specific cases where 
solutions can be found analytically. There are also some 
general qualitative features of the dynamics of morpho- 
logical evolution which can be deduced from the form 
of equations (|J). The right-hand side of equation (2b) 
has three negative and one positive terms. The latter, 
which is twice the expected number \xh of new traits per 
unit time interval per lineage, gives the maximum possi- 
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ble rate of increase in morphological disparity D. Each 
of the three negative terms is proportional to D mean- 
ing they are negligible initially when the clade is confined 
to a small volume on the hypercube (when D is small) 
but become increasingly important as the clade diversi- 
fies morphologically (when D becomes larger). The first 
negative term in the right-hand side of equation (2b) is 
related to the geometric structure of the morphospace: 
as the clade expands in the morphospace, it becomes less 
and less probable that a random morphological change 
will lead outside the volume of the morphospace already 
occupied by the clade. The second term describes the re- 
duction in D because of the splitting of lineages into in- 
dependent units without immediate and significant mor- 
phological changes. The third term specifies the reduc- 
tion in D due to the extinction of subclades. Thus, equa- 
tion (2b) predicts rapid initial increase in D with a slow- 
down afterwards coming because of the geometric struc- 
ture of the morphospace and the effects of extinction and 
speciation on morphological disparity. The slowdown is 
expected even when all underlying processes are time- 
homogeneous. 

Let us turn now to variable d. The general solution of 
equation (2c) can be approximated as 

d = - ( 1 - exp(- J 2fidt) \ . (3) 

This shows that provided /i > 0, the average distance d 
of the clade from its ancestral state is always expected 
to increase monotonically towards the asymptotic value 
L/2. In particular, this dynamical feature is not affected 
by changes in the rates of extinction, origination and mor- 
phological changes. 

To illustrate the dynamics of the clade diversification 
let us assume that all parameters of the model are con- 
stant and that the rate of origination exceeds that of ex- 
tinction. Let R = a — 8 be the clade's growth rate. In 
this case, the clade size increases exponentially: 

N(t) = exp(r), (4a) 

the average pairwise distance within the clade tends to 
L/(2 + 5 2 ((>/2fj i y. 

D(t) = Le- aT exp(be- T )cb a [F(-a,6e" T ) - T(-a,b)] , 

(4b) 

whereas the average distance from the founder approaches 
L/2: 

d{t) = - [1 - exp(-CT)} . (4c) 

Here r = Rt,a= (4/i + S 2 (f)/R,b = 2a 1 /R,c = 2fx/R 
and r(x, y) is the incomplete gamma-function (e.g. Grad- 
shteyn & Ryzhik, 1994). The dynamics effectively depend 
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Figure 1: Dynamics of diversification as predicted by equa- 
tions (4) with a = 0.32,5 = 0.25,^ = 0.006 and the unit 
time interval corresponding to one million years (cf. Foote 
1996a). (a) Zero probability of speciation of the first type 
(<7i = 0, (72 = 0.32), (b) Equal probabilities of speciation of the 
first and second type (<n = 0.16, 02 = 0.16), (c) Zero probabil- 
ity of speciation of the second type (ci = 0.32, 02 = 0). Thick 
lines represent clade size N relative to the size achieved by the 
end of the time interval studied, N max = 1096. Dashed lines 
represent the average distance of the clade from the founder 
relative to the asymptotic value L/2. Thin lines show dis- 
parity D (normalized relative to L/2) corresponding to the 
extinction every 10-th time step of a sub-clade representing 
10% (top lines), 20% (middle lines) and 30% (lower lines) of 
the clade. 

on only three parameters, (J./R, 52<p/ R and <7i/R, charac- 
terizing the probabilities of a morphological change, sub- 
clade extinction and speciation of the first type relative to 
the clade's growth rate. Figure 1 illustrates the patterns 
of diversification predicted by equations (Q) . The numer- 
ical values used for the overall extinction and origination 
rates 6 and a are the same as in Foote (1996a) where the 
estimate of S was based on Raup's (1991) data whereas 
the origination rate a was set to produce an increase in 
diversity to about 1000 lineages in 100 million years. The 
numerical value used for \i was estimated from the blas- 
tozoan data (see below). Depending on parameter value, 
disparity D can increase faster or slower than d. Note that 
even in the case of exponential increase in the clade size, 
the clade's disparity D will approach L/2 asymptotically 
only if there is no sub-clade extinction (62 = 0). In other 
situations, the asymptotic value of D will be smaller than 
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L/2. For instance, if on average every 50-th time step a 
subclade representing 50% of the clade goes extinct, then 
5 2 = 1/50x1/2 = .01,0 = 1/2, and if /a = 0.0025, then D 
will approach L/3. The Appendix lists several other spe- 
cific cases where equations (2) can be solved analytically. 
In general, if the clade increases in size (if R > 0), its mor- 
phological diversification is the fastest initially and slows 
down afterwards. The dynamics of D are close to that 
in the exponential case, especially during initial stages. 
If the clade decreases in size (if R < 0), its loss of mor- 
phological disparity is delayed relative to the loss in the 
number of lineages (cf. Foote 1993, 1996a). 

3. APPLICATION OF MODEL 

The model presented above provides a framework for 
studying the complex processes of clade diversification. It 
can be used to train our intuition about these processes, 
to identify key components, and to suggest hypotheses 
that can be tested against fossil data. Can one use the 
observed dynamics of D and d for making quantitative in- 
ferences about the underlying processes? Using disparity 
D is complicated because both its dynamics and asymp- 
totic value depend on a number of parameters that may 
change in time and be difficult to extract from fossil data, 
such as the relative importance of clade - vs. species - ex- 
tinction, the relative importance of phyletic vs. cladoge- 
netic change, and the proportion of speciation events that 
involve significant morphological change. In a wide vari- 
ety of curcumstances, however, distance d is expected to 
increase regularly towards a fixed asypmtotic value L/2. 
Equation (§) implies that -ln(l - 2d/L) = J^fxdt. If 
the overall probability of a morphological change /i does 
not change (significantly) in time, the integral is simply 
2/it and the dependence of q = 1 — 2d/L on time should 
be a linear function on the semilog-scale. (In statistical 
physics variable q is known as the average "overlap" be- 
tween binary sequences, e.g. Derrida and Peliti 1991.) 
This provides a simple test for approximate constancy 
of (i. The constancy of \x during a certain period would 
suggest the constancy of evolutionary rates and mecha- 
nisms. Thus, if the rate of increase of disparity D declines 
over time while the rate of increase of —ln(q) remains ap- 
proximately constant, then the explanation of the decel- 
eration of morphological diversification as a consequence 
of a change in evolutionary rates should be rejected. If 
— ln(q) changes as a linear function of time, the slope of 
the regression line gives an estimate for 2/i. In our model, 
the overall probability of a morphological change /i is a 
sum of the phyletic component (i± and the cladogenetic 
component [12^2- Given some information about the pat- 
terns of origination in a clade, it should be possible to help 
constrain the relative importance of phyletic and cladoge- 
netic change. For example, a reduction in 02 is supposed 



to translate into a comparable reduction in /i unless the 
phyletic component is much larger than the cladogenetic 
component. Thus, if there is evidence for a change in spe- 
ciation rate (a) without a proportional change in the rate 
of morphological evolution (/i), this suggests a greater 
role for phyletic evolution than would be suggested by a 
concordant changes in fi and a. The methods of extract- 
ing various information about /1 from the clade level data 
proposed above are potentially useful in general, but es- 
pecially in those cases in which the phylogenetic informa- 
tion needed to measure ancestor-descendant differences is 
unavailable. 

4. DIVERSIFICATION OF BLASTO- 
ZOANS 

In light of the foregoing discussion, I use the model de- 
scribed above to reanalyze the data on the morphological 
diversification of blastozoans (Foote 1992, 1996a). These 
data provide one of the best illustrations of the pattern 
of accelerated early morphological diversification. The 
data represent sixty-five discrete characters measured for 
147 species spanning across 12 stratigraphic levels from 
the Lower Cambrian to the Permian. Morphological dis- 
tances between species were measured as the total number 
of differences divided by the number of characters com- 
pared corresponding to d/L and D/L in the notation of 
the model. The time scale used is from Tucker & McK- 
errow (1995) (with the exception of the Carboniferous 
and Permian which are not covered by Tucker & McKer- 
row and for which I used the Harland et al. 1989 scale). 
The data points were placed in the middle of the inter- 
vals. Morphological disparity (D) increases more rapidly 
than taxonomic diversity (N) reaching one half of the 
maximum observed level by the Late Cambrian and the 
maximum observed level by the Middle Ordovician (Fig- 
ure 2a, b). As emphasized earlier (Foote 1992, 1996a), 
during the period of this clade's expansion the rate of 
increase of D apparently declines in time (Fig. 2b). The 
taxonomic diversity N grows through the Caradocian and 
decays after it suggesting a major change in the pattern of 
origination and extinction somewhere near the Caradoc- 
Ashgill boundary. This change coinsides with apparent 
drops in both d and D. A drop in morphological dis- 
parity D can be caused by an increase in subclade ex- 
tinction rates, an increase in the rate of speciation of 
the first type, the selective extinction of morphologically 
"peripheral" (relative to the founder) lineages, or the se- 
lective proliferation of morhologically "central" lineages, 
among other factors. A drop in d can be caused by extinc- 
tion of morphologically "peripheral" lineages and/or by 
intensive speciation of "central" lineages. The overall de- 
crease in taxonomic diversity N between the Caradocian 
and Ashgillian suggests the possibility that it may have 
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been increased extinction of morphologically peripheral 
lineages that caused the drop. Note that the decrease in 
morphological disparity is delayed relative to the decrease 
in taxonomic diversity. Distance d continues to increase 
after the Ordovician while D declines. The fact that D 
is low and d is high later in the clade history means that 
the clade forms a compact group evolving far away from 
the founder. (In this case it is the Blastoidea [Foote, 
personal communication].) This does not necessarily im- 
ply directionality in the processes governing clade evolu- 
tion or species selection. Such behavior is expected for a 
completely random walk in a multidimensional space (cf. 
Charlesworth 1984; Bookstein 1987). 

As indicated by the apparent linearity on the semilog- 
scale (Figure 2c), the dynamics of d appear to be time- 
homogeneous from the Lower Cambrian through the Mid- 
dle Ordovician and from the Upper Ordovician through 
the Upper Carboniferous with the drops in d near the 
Middle Ordovician-Upper Ordovician boundary and in 
the Permian. I used this as a justification for splitting 
the data set into two parts (for computing separate re- 
gression lines) and excluding the Permian point. The 
estimates of /i for the periods from the Lower Cambrian 
through the Middle Ordovician and from the Upper Or- 
dovician to Upper Carboniferous are (5.8 =F 0.4) x 10~ 3 
and (3.6 =F 0.2) x 10 -3 , respectively. For exponential pro- 
cesses, the time scale is usually characterized in terms of 
a half-life T 1/2 . The half-life for d is ln(2)/(2/x). With 
/i = 0.0058 and fi = 0.0036, T 1/2 is about 60My and 
96My, respectively. The linear regressions provide an ex- 
celent fit (the r 2 coefficients are 0.983 and 0.986), the 
slopes are significantly different from zero (P < 0.001) 
and from each other (P < 0.01). The increase in the 
quality of fit gained by fitting two separate lines rather 
than just one line is significant at P < .01 (the Sncdecor 
test). The numerical values reported above should be 
taken with some degree of caution for some of the as- 
sumptions underlying the regression methods might be 
violated. The regression estimates are obviously sensitive 
to details of time scale and resolution. The appearance 
of more precise stratigraphic data would probably require 
to reevaluate parameter estimates but is not expected to 
change our qualitative conclusions about the diversifica- 
tion of blastozoans. 

Overall, the data are compatible with a moderate 
(38%) reduction in fi that took place near the bound- 
ary of the Middle Ordovician and the Upper Ordovician. 
A decrease in the rate of morphological evolution was also 
advocated on the basis of the shape of the disparity curve 
(Foote 1992) and on estimates of morphological separa- 
tion between closely related taxa (Wagner 1995a). The 
reduction in fi, however, is not responsible for the appar- 
ent deceleration of morphological diversification observed 
during the first third of the clade's history when \x was 
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Figure 2: Reanalysis of the data (Foote 1992, 1996a) on the 
morphological and taxonomic diversification of Blastozoans. 
(a) The number of genera (which is assumed to be a reason- 
able proxy of the number of species N). (b) Morphological 
disparity and the average distance from the founder. For each 
stratigraphic level except the first one, the average distance 
from the species-founder d was approximated by the average 
distance between the current level and the first level. For the 
first stratigraphic level, d was approximated as half of D. (c) 
Transformed values of d and the corresponding linear regres- 
sion (see text for more details). In Figure (a) the 14 data 
points correspond to the following stratigraphic levels: Cam- 
brian (Lower, Middle/Upper), Ordovician (Tremadoc, Arenig, 
Llanvim, Llandeilo, Caradoc, Ashgill), Silurian, Devonian 
(Lower, Middle/Upper), Carboniferous (Lower, Upper) and 
Permian. Figures (b) and (c) use a coarser resolution with 
only 3 Ordovician intervals (Lower, Middle and Upper). 

apparently constant. 

The fossil record shows that there has been significant 
decline in rates of origination within major taxa through 
their histories (e.g. Van Valen 1985; Sepkoski 1998). The 
moderate size of the decrease in /j, observed for blasto- 
zoans together with a significant decrease in speciation 
rates would suggest that for this clade morphological evo- 
lution is driven mainly by anagcnetic rather than clado- 
genetic changes. Although the fact that blastozoan tax- 
onomic diversity increased initially and declined later on 
(see Fig. 2a) is compatible with the decline in speciation 
rates this is not a definite conclusion. The decline in tax- 
onomic diversity can be caused by an increase in extinc- 
tion rates rather than by a decrease in origination rates. 
Additional data are needed for reaching more precise con- 
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elusions. 

5. DISCUSSION 

It is important to realize that apparent secular changes in 
the rates of morphological evolution for a clade as a whole 
do not necessarily mean secular changes in the processes 
acting at the level of individual lineages. In particular, 
the observed deceleration of morphological disparity does 
not necessarily imply a decline in the size or probabil- 
ity of morphological changes for individual lineages. The 
model presented here has demonstrated that such a de- 
celeration is expected from the geometric structure of the 
morphospace and the effects of extinction and speciation 
on morphological disparity even when all relevant pro- 
cesses are time-homogeneous. These theoretical predic- 
tions appear to be very robust. In particular, the differ- 
ences between exponential and logistic growth in the tax- 
onomic diversity do not translate into significant changes 
in the corresponding dynamics of morphological diversi- 
fication. Our basic conclusions will definitely be valid if 
the traits have more than two discrete states and should 
be valid if the morphological space is continuous rather 
than discrete as long as it is finite. In the continuous 
case, the effects of speciation of the first type and sub- 
clade extinction on disparity will be similar to that in the 
discrete model considered here. Although demonstrating 
the existence, nature, and importance of morphological 
boundaries can be difficult (McShea 1994, Foote 1996a), 
it is intuitively obvious that, given they exist, the pro- 
cess of divergence will slow down even if these two fac- 
tors are absent. The potential importance of the model 
I present is not only that it quantifies and trains our in- 
tuitions, but also that it allows one to test whether, in 
the case of discrete characters, the observed deceleration 
in morphological diversification is likely to be a simple 
consequence of the nature of evolution on a binary hyper- 
cube, and thus whether it is necessary to invoke tempo- 
ral heterogeneities in evolutionary rates and mechanisms 
to explain an observed pattern. The model makes fal- 
sifiable predictions about the dynamics of morphological 
disparity and the average distance of the clade from its 
ancestral state, provides a simple method to evaluate the 
rate of morphological evolution and suggests an approach 
for comparing the importance of anagenetic and cladoge- 
netic changes in morphological diversification. In the case 
of blastozoans, I find no evidence that major changes in 
evolutionary rates and mechanisms are responsible for the 
deceleration of morphological diversification seen during 
the period of this clade's expansion. At the same time, 
there is evidence for a moderate decline in overall rates 
of morphological diversification concordant with a major 
change (from positive to negative values) in the clade's 
growth rate. 



The model has its limitations. The most significant 
is probably that it describes only average behavior and 
says nothing about variation which will always be present 
in the fossil record (and numerical simulations). In par- 
ticular, this makes it difficult to evaluate the statistical 
power of the test of time-homogeneity proposed above. 
The model developed above makes no restrictions on mor- 
phology in the sense that all character combinations are 
assumed to be potentially realizable. In terms of the 
metaphor of "adaptive landscapes" (Wright 1932), the 
model assumes a "flat" landscape similar to those in mod- 
els of neutral molecular evolution (e.g. Derrida & Peliti 
1991). In general, because of genetic, developmental, 
or ecological constraints some of the possible character 
combinations can be prohibited. In this case, the mor- 
phospace will be mathematically equivalent to a hyper- 
cube with "holes" (with "holes" representing prohibited 
character combinations) and the corresponding adaptive 
landscape will be "holey" (Gavrilets 1997, 1999; Gavrilets 
& Gravner 1997; Gavrilets et al. 1998) rather than "flat" . 
If the proportion of holes is not extremely high, "viable" 
character combinations will form a "giant" cluster extend- 
ing through the whole morphospace. A characteristic sig- 
nature of a random walk on the giant cluster appears to be 
a stretched exponential dependence of overlap q on time 
(e.g. Lemke & Campbell 1996): q(t) ~ ea;p(-(t/T) /3 ), 
where r and (3 < 1 are parameters (with no holes f3 = 1). 
The fitting of the stretched exponential curve to blasto- 
zoan data has led to inconclusive results: although the fit 
is good, it is not better than the fit of a simple exponen- 
tial curve described above. More detailed data sets are 
needed for more precise conclusions. 
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Appendix 

Each new trait in a lineage increases or decreases its dis- 
tance from the founder by one with probability 1 — d/L 
and d/L, respectively. Extinction of a lineage is not ex- 
pected to change d. Thus, the overall expected change 
in d is u[l x (1 — d/L) — 1 x d/L] which reduces to (2c). 
Consider the dynamics of D. Because each lineage in 
a randomly chosen pair can evolve morphologically, the 
rate of change in D induced by phyletic evolution and 
by cladogenesis is twice as big as in the case of d and is 
—4fi(D — L/2). Speciation events of the first type will 
decrease D because the distance between the immediate 
descendants of a species that has split will be zero. Let 
P be the probability that two randomly chosen lineages 
originated from a split of a species in the previous time 
interval. If k and var(k) are the average and the variance 
of the number of "offspring" species that a species leaves 
(counting itself) in the next time interval, then probabil- 
ity P can be represented as P « (var(k)/k + k — 1)/N 
(Crow & Kimura 1970). Assuming that a proportion S of 
species go extinct whereas the remaining species survive 
and speciate, k = 1 + a — 5, var(k) = a + S where the 
last equality assumes that both extinction and speciation 
rates are small. Thus, P « 2a/N (cf. Derrida & Peliti 
1991). Distance D after the speciation events of the first 
type can be represented as (1 — Po\/o) x D + Pa\/a x 0. 
Thus, the expected reduction in D due to the speciation 
events of the first type is —2a\/N. Extinction of individ- 
ual lineages is not expected to change D. Let us assume 
that there are / T-subclades and let Df, be the average 
morphological distance between two T-subclades. The 
average distance within the whole clade can be approxi- 
mated as (N/f) 2 Dbf(f — l)/N 2 . Accordingly, the average 
distance after extinction of a T-subclade is approximately 
(N/f) 2 D b (f - 1)(/ - 2)[(1 - l/f)N} 2 . This shows that 
extinction of a T-subclade reduces D by D/ f 2 . Let n be 
the probability of extinction of a T-subclade per unit time 
interval and (j> — 1// be the proportion of the clade that 
goes extinct. Then the rate of reduction in D due to sub- 
clade extinction is approximately —52<j) where 82 = rj<f>. 
Thus, the overall change in D is given by equation (2b). 

In the main text, a case with all parameters constant 
was considered. Here I list several other cases that can 
be treated analytically. Let the growth rate decrease lin- 
early with the clade size: R = a — 8 = r(l — N/K). Then 
the clade size approaches the "carrying capacity" K ac- 
cording to the logistic curve N(t) = Ke T /(e T + K — 1) 
where t — rt. The difference between exponential and 
logistic growth in N should be most important after a 
transient time when the exponential model predicts very 



large values of N whereas in the logistic model N ap- 
proaches the carrying capacity K. However if N is large, 
the second term in the right-hand side of equation (2b) is 
negligible. Thus, the difference between exponential and 
logistic models for N is not expected to translate into 
significant changes in the dynamics of morphological evo- 
lution if K is not too small. Below I make this argument 
more precise. The decrease in the growth rate R with 
the clade size can result from decrease in the origination 
rates and/or increase in the extinction rates. I assume 
that other parameters do not change. If 02 = const, then 
the dynamics of d are still described by equation (4c). 
If o~2 decreases linearly with the clade size N from 172(1) 
to a%{K), then the dynamics of d are approximated by 
equation 



-,-Ct 



e 

1+ K 



(5) 



Here C = 2fi*/r,fi* = Mi +^2(1), Aa 2 = a 2 (l) -a 2 (K) 
is the overall change in a 2 , and it is assumed that K » 
1. To solve equation (2b) analytically when extinction 
and/or origination rates change with N one needs ad- 
ditional simplifying assumptions. If 5\ changes with N 
whereas all other rates are constant, the dynamics of D 
are described by (4b) with a = (4[i + 2cri/K + 52(f))/r, b = 
2(1 — X/K)ai/r, c = 2/i/r. Note that if the "carrying ca- 
pacity" K is large, coefficients a, b and c are close to the 
values corresponding to the exponential growth case and 
the dynamics of D under exponential and logistic growth 
are similar. If o\ decreases linearly with the clade size 
N whereas other rates are constant, the D evolves ac- 
cording to (3b) with a = (4/i + 2(cn(l) - r)/K)/r),b = 
2(1 — 1/ K)cri(l) /r, c — 2/x/r. If the growth rate r is small 
relative to the origination rate ci(l), the dynamics of D 
will be similar to that under exponential growth. If all 
speciation events are of the second type (<ti = 0) and 
there is no family extinction (<5 2 = 0), then the dynam- 
ics of D are described by the right-hand side of equation 
(||) with C = 4/x*/r. Let the clade size decrease linearly 
in time: N(t) = N(0) — rt. Then the average pairwise 
distance D changes according to equation 

+ca b - 1 e aN{t) N{t) b [T(l - b, aN{t)) - F(l - b, aN(0)} 

where r = rt, parameters a, b and c are defined below 
equation (4b), and D(0) is the initial value of D. 



